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Magnetic Moments of A and Q Baryons with Dynamical Clover Fermions 
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We calculate the magnetic dipole moment of the A(1232) and f2~ baryons with 2+1-flavors of 
clover fermions on anisotropic lattices using a background magnetic field. This is the first dynamical 
calculation of these magnetic moments using a background field technique. The calculation for 
is done at the physical strange quark mass, with the result in units of the physical nuclear magneton 
/^Q- = — i.93(8)(12) (where the first error is statistical and the second is systematic) compared 
to the experimental immber: -2.02(5). The A has been studied at three unphysical quark masses, 
corresponding to pion mass = 366, 438, and 548 MeV. The pion mass dependence is compared 
with the behavior obtained from chiral effective field theory. 

PACS numbers: 11.15.Ha,12.38.Gc,13.40.Em,14.20.Gk 



Calculations of hadron properties from first principles 
using lattice QCD have been rapidly advancing in re- 
cent years. The newly available fully-dynamical (un- 
quenched) lattice configurations have made it possible 
to significantly reduce the systematic error of lattice cal- 
culations. Of the properties that can now be reliably 
computed on the lattice are the electromagnetic (e.m.) 
properties of baryons and in particular their electromag- 
netic moments. Here we present a first dynamical cal- 
culation of the magnetic dipole moment of the A(1232) 
and f2~(1672) baryons using a background cm. field. 

These particular baryons are chosen for the following 
reasons. They both are distinguished members of the 
baryon decuplet and as such they have much in common. 
On the other hand, while the magnetic moment of the 
is measured to a few-percent accuracy, the tiny lifetime 
of the A resonance (« 6 x 10"^'' sec) hinders the deter- 
mination of its magnetic moment and the experimental 
efforts are still ongoing. This is why a simultaneous lat- 
tice calculation for these two baryons can both be tested 
against experiment in the case and provide predictions 
in the case of the A. 

The background-field method adopted here is presently 
the simplest and cleanest way to access the static e.m. 
moments on the lattice, as it amounts simply to mea- 
suring the shift in the mass spectrum upon applying 
a classical background field [l| (for most recent appli- 
cations to baryons in the quenched approximation see 
Rcfs. The other possibility is the form factor 

calculation extrapolated to the — point (from the 
minimum momentum-transfer on the lattice, which is 
27r/aL, with L the number of points in the spatial direc- 
tion). However, in comparison with the background- field 
method, this method is additionally complicated by the 
noise of the three-point function calculation as well as the 
uncertainties in the extrapolation, see Refs. [1, 0] for 
recent calculations of the A e.m. form factors (the first 
calculation is done in the quenched case, and the latter 
in the dynamical case). 

In order to calculate the magnetic dipole moments, we 



implement the constant background magnetic field in the 
following fashion. On a given configuration, we multiply 
all of the S'J7(3) gauge fields by a U{1) gauge field, and 
invert the Dirac operator on that background to get the 
quark propagator in a background field. The f/(l) links 
are given by 



Ufiix) = exp [iqaA^{x)] 



(1) 



where q is the charge of the quark whose propagator we 
are calculating. For a constant magnetic field with a 
magnitude of B pointing in the -l-z-direction, the usual 
choice is A^{x, y, z, t) = aBx 6^y. The problem with this 
choice is that due to the condition that the gauge links 
must be periodic, the field is continuous only if qa^B = 
iim/ with integer n. Hence, the minimal value of B is 
severely limited by the size of the lattice. This limitation 
is somewhat relaxed for the following choice of the field ml, 
0J: 



A^{x,y,z,t) 



aBx 5^y 
-aBLy 5^^ 



iix^L-1 
iix^L-1 



(2) 



Thus, all of the y-links are modified by exp \^^qa?Bx^^ , 
all x-links on the x boundary are modified by 
exp [—iqa^BLy] , and all other links are unchanged. The 
additional modification of the links on the a;-boundary, 
allows us to achieve continuous constant field everywhere 
on the lattice with a more relaxed constraint on the value 
of the field: 



qa'^B 



IF' 



(3) 



The latter periodicity constraint corresponds with the 
more physical requirement that the magnetic flux (pla- 
quettc) remain continuous through the boundary. 

With this 1], one can calculate a baryon two-point 
function which behaves for large time in the usual manner 



C(t) ~Ae-™(-^)*-h... 



(4) 
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but with the exponential damping governed by a _B-field 
dependent mass 



m{B) =mB- ^^zB + 0{B^ 



(5) 



Here, ms is the mass of the baryon in the absence of any 
external field. The magnetic moment along the direction 
the field is given by 



^J■z = fJ-Sz/S , 



(6) 



where ^ is the value of the magnetic moment, Sz is the 
spin projection and S the total spin (in our case — 3/2 
and so Sz can take the values ±3/2 and ±1/2). It is 
useful to form the quantity 



Ams, = m{-B) - m{B) = 2fiB^ + 0{B^) 



(7) 



to cancel the effect of the next order in the B expansion. 
We have computed this quantity for all the various spin- 
projection values and extracted the magnetic moment 
from the following combination: 
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8S 



A' 



"^3/2- 



' Am 



-3/2 ^ 



3(Ami/2- Am_i/2) 



(8) 



Since the Amg^ are highly correlated, the error on /i 
is determined using a jackknife, and this combination is 
chosen to average the results from all spin components 
under the jackknife procedure. 

On a technical note, the number input into the simula- 
tion is qa^B, and thus includes the product of the quark 
charge and the magnetic field in lattice units. In order 
to account for the quark charges of the up, down, and 
strange quarks, for a single magnetic field B we must 
use two values of qa^B, corresponding to the fact that 
Qu = —2qd,s- For particles made up with only a single 
quark (A++,r2~), we need one input value, but for the 
A"*" or the nucleon, for example, we must use two in- 
puts that differ by a factor of —2 so that the quarks all 
experience the same S-field. 

Notice that even with the modified periodicity con- 
straint in Eq. ([3]), the minimum value of the magnetic 
field may still be large enough to distort the baryons, 
and thus introduce errors into the extrapolation of the 
magnetic moment. Ideally we would use volumes large 
enough that this would not be true, however this can be- 
come rather expensive. Earlier studies [l|, 0] , ignore the 
periodicity constraint,^ using small fields that would not 
distort the particles and also ensure the linear relation- 
ship between the extracted mass in Eq. ([5]) and the mag- 
netic field. In addition, they imposed Dirichlet boundary 



^ These studies do not include the modification of the x-Unks on 
the boundary, and thus their periodicity constraint is given by 
qa?B = 2nn/L, an order of magnitude larger than ours. 



TABLE I: Lattice parameters used for the current work. 
On these lattices, the anisotropy is at /as ~ 3.5, and we 
show the pion masses on the lattices in physical units. The 
<ia?7ivai ~ —0.0743 dataset was used to calculate the Q~ mag- 
netic moment, so the mA,n listed in that row is the fl~ mass, 
while the other baryon masses are those of the A. 



Volume 




nin 


mA,n # 






(GeV) 


(MeV) (GeV) configs 


16^ 


X 128 


-0.0808 


548 


1.562 110 


16^ 


X 128 


-0.0830 


438 


1.485 91 


24^ 


X 128 


-0.0840 


366 


1.408 202 


24^ 


X 128 


-0.0743 


366 


1.65 213 



conditions and place the source in the center of the lat- 
tice to perhaps ensure the quarks will not feel the effects 
of the discontinuity. 

The difficulty with this approach, however, is that 
there are significant finite volume effects in the results 
in the magnetic moments, using their implementation. 
Specifically, in the quenched calculation of Ref. the 
authors see effects that are as large as 35% for the light- 
est pion mass when comparing the 16'^ volume to a 24'^ 
volume, with a lattice cutoff of a^^ w 2 GeV, and a 
pion mass of about 522 MeV. Since taking the pion mass 
closer to the physical point, finite volume errors become 
more substantial, we would like to reduce the finite vol- 
ume effects coming from the background field as much as 
possible. 

Using the implementation of the magnetic field above, 
we have shown that ignoring the periodicity constraint 
somewhat will not introduce noticeable finite volume er- 
rors coming from the background field, so long as one 
uses the implementation of the background field shown 
in Eq. ^ p\. With this method, we are able to trust 
results coming from simulations on smaller volumes (still 
keeping Tti-j^L ^ 4 so wG can minimize finite volume effects 
coming from a small pion mass), which are less expen- 
sive. Using the methods in Ref. 3], for example, one is 
restricted only to larger volumes. 

We now present our results, which are the first dy- 
namical calculations for magnetic moments using a back- 
ground field. 

We use dynamical anisotropic lattices with 2+1 flavors 
of Stout-smeared Clover fermions [^, [l^ , on two volumes 
and a single lattice spacing. We show the relevant param- 
eters in TablelD Note that on these lattices, a bare quark 
mass parameter of —0.0743 corresponds to the physical 
strange quark mass. Both sets of configurations have 
an inverse spatial lattice spacing of 1.61 GeV and an 
anisotropy of about 3.5 (so a^^ w 5.61 GeV). More com- 
plete information on these configurations, specifically the 
tuning of the lattice parameters, can be found in [3, [lot . 

In all A cases we used three magnetic fields for the 
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simulations, corresponding to n = ±1/2, ±1, and ±2 in 
Eq. ([3]). For the fi, we only used n = ±1/2 and n = 
±1. The n = ±1/2 field does not satisfy the periodicity 
constraint. We expect the errors entering here due to 
this to be negligible as was shown in Ref. [8|] . Even with 
these fields, we already see higher-order terms appearing 
in the expression for the masses extracted from two-point 
functions, so larger magnetic fields will begin to introduce 
effects coming from even higher-order terms in Eq. ([8|. 
We calculate all four spin projections for the baryons, as 
well as using both positive and negative magnetic fields, 
and we average over all of these to reduce the errors. 
Additionally, on each configuration, we calculated the 
quark propagators starting from four time sources t = 
0, 32, 64 and 96, using the EigCG algorithm developed 
in Ref. ll| to decrease significantly the time it takes to 
invert the Dirac operator. 

In Fig. [1] we show the A+^" magnetic moments in units 
of the physical nuclear magneton fiN, for the three pion 
masses simulated. One can see noticeable effects coming 
in at 0{B'^), and we have fit each dataset to a quadratic 
form 



4.5 



/i = /io + b{ea^B) 



(9) 



With each set, since we only have two parameters and 
three data points, a correlated fit is not possible, and we 
do not take the fitted value of (in principle this would 
be the most appropriate value, as it subtracts out the 
dependence) and its error as our final result. Instead, we 
see that for the smallest of the i3-fields simulated, the 
0{B^) effects are small, and so we take that data point 
as our determination of the magnetic moment (in fact, 
the data point at the smallest value of B is consistent, 
as one can see from the Figure, with the value of /ig). 
The fits performed give drastically smaller errors than 
the data, and so we choose to use the error from the 
data, to account for possible uncertainties in this method. 
There is a slight shift (within errors) in the extracted /iq 
compared with the smallest i3- field data point, and for 
the heavier mass this is its largest at 5%. Similar results 
can be seen for the A+ and 17 ~, and we show all of our 
results in Table. HIl Note that our results for the A~ are 
not included, because with this method of calculation, 
we have the exact equality /iA++ = ~2/iA-. 

Also in the table, we show the experimental numbers 
for these quantities. We can see that for the A"*""*", there 
is a possible upward trend as we decrease the pion mass. 
We illustrate this in Fig. [2l where we show the present 
results together with the chiral effective field theory cal- 
culations of Ref. [3] for the dependence of /^a+- 
These calculations have one free parameter for /iA+, cor- 
responding with its value in the chiral limit. We also 
indicate a theoretical error band, corresponding with an 



next chiral order, with m^y^ the physical pion mass and 



error of (m7r+mP^'')/m^"'"' estimating the corrections of 
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FIG. 1: In this figure we show the magnetic moment in units 
of the physical nuclear magneton for the three input magnetic 
fields used, as well eis quadratic fits to each dataset, to remove 
residual dependence in the magnetic moments. 



TABLE 11: Calculated magnetic moments in units of /ijv, the 
physical nuclear magneton (taken as the value for the data for 
the smallest B-field, as discussed in the text). For compari- 
son, we have combined all experimental errors in quadrature. 





Ma++ 




Mao 


Mn- 


548 


3.65(13) 


2.60(8) 


-0.07(2) 




438 


3.55(14) 


2.40(5) 


0.02(3) 




366 


3.70(12) 


2.40(6) 


0.001(16) 


-1.93(8) 


PDG: 


5.6(1.9) 


2.7(3.5) 




-2.02(5) 



m 



phys 
A 



the physical A mass. One notices a strong cusp 



behavior for the real part of /ia , which is due to the open- 
ing of the A — > TriV decay channel. Therefore, no strong 
conclusions on the value of /ia at the physical point can 
be made until this extrapolation has been done. We leave 
such a systematic study for a future work. On the ex- 
perimental side, there are new experiments from MAMI 
for the magnetic moment of the A+ with much reduced 
errors, yet these have not yet been fully analyzed. 

Since the sea quarks do not carry electric charge (which 
is the case for all current lattice simulations), there is a 
relationship that holds within the quark model in the 
isospin limit, where /iA++ — 2/iA+- Clearly this relation- 
ship does not hold with our results above, but we could 
use this relationship to reduce the systematic uncertain- 
ties in our determination. This would clearly increase the 
values obtained for the A++ and reduce it for the A+. 

As for the fi^, the strange quark mass is close to its 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 
[GeV^] 

FIG. 2: Chiral effective field theory calculations of Ref. [l^ 
for the m,r dependence of , in units of the physical nu- 
clear magneton. Both real and imaginary parts of /^a are 
displayed. For the former, the bands show a theoretical error 
estimate, as described in the text. The value at the physical 
pion mass corresponds with the experiment of Ref. tl3| , where 
both statistical and systematic errors are displayed. 



physical value (as we can see by the fact that the Vt~ mass 
is close to the observed value), so we expect the result 
to match more closely to the experimental value. As we 
can see, it agrees tremendously well. This agreement 
is expected, as quantities involving the should have 
little dependence on the light sea quark mass. On the 
two B-fields we simulated for the 0~, we see a slight 
dependence in the magnetic moment, roughly of the same 
size as for the A. Additionally, we see that the errors 
associated with the experimental value are comparable 
to the statistical lattice errors here. 

In order to improve on the quoted results, we must ac- 
count for the systematic errors that arise from a variety 
of sources in the calculation. First there is the finite lat- 
tice spacing, which is difficult to estimate given the lack 
of any calculations of the magnetic moments (quenched 
or dynamical) at multiple lattice spacings. Given we are 
using Clover fermions, errors of 0{a) disappear, so one 
would expect errors to be roughly 0[a?K^^Y)) ^0.03. As 
it will be some time before a second lattice spacing is 
available on these configurations, we will assume there is 
a 3% systematic error that arises from the finite spatial 
lattice spacing. 

Additionally there are errors arising from remnant fi- 
nite volume effects, coming not from the background 
field, but from the pion mass. These are most likely 
negligible since in all cases, m^ri ^4, and thus the errors 
from these effects are less than e~™"^ « 2%. 

Finally, there are uncertainties that plague any current 



calculation of the magnetic moments on the lattice, being 
that the sea quark charges are set to zero. We expect 
these to not be very large, coming from the discarded 
diagrams in which the valence quarks in the baryons emit 
photons that couple to the sea quarks. These at most are 
of order a relative to the terms that are included, and 
thus are expected to be at most 1%. Related errors are 
those coming from the extrapolation, and this is going 
to be at most 5%, as mentioned below Eq. 

We have presented here the first (using a background 
field method) dynamical results for the A and fJ" mag- 
netic moments on dynamical 2+1-flavor lattices, which 
are consistent (given the pion mass used) with experi- 
mental values that have been measured. Presently, the 
accuracy obtained in the lattice result for the f2~ mag- 
netic dipole moment is comparable with the experimental 
accuracy. We can use the above discussion to estimate 
the systematic error on our result for the magnetic 
moment. We make a conservative estimate, and use the 
maximum values for each source of systematic error, and 
add those in quadrature, giving an error of 6%. Thus we 
quote jjL^- — — 1.93(8)(12)^7v for our final result. 

To make significant progress on these results, simula- 
tions going to lighter pion masses, especially below the 
A — > t:N threshold, are essential in precisely determining 
the magnetic moments. Nevertheless, it is rather encour- 
aging that one can obtain already such precise results 
with the resources currently available. 
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